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1 .  INTRODUCTION 


It  Is  vsll  known  chat  any  raal  continuous-data  function  can  be  represented 
by  a  Fourier  series  of  infinite  tens,  provided  a  certain  set  of  conditions 
are  net.  In  practice,  the  infinite  series  is  truncated  to  contain  only  a 
finite  number  of  terms.  Better  approximation  is  obtained  if  more  terms  are  in¬ 
cluded  in  the  series.  This  last  statement  is  not  exactly  true  for  the  case  of 
a  real  discrete-data  function.  For  this  case,  there  is  an  optimum  truncation 
for  its  Fourier  series  representation.  This  fact  has  not  been  widely  recog¬ 
nised  by  practicing  engineers.  This  report  discusses  this  Important  and  inter¬ 
esting  fact. 

II.  DISCUSSION 

A  real  continuous-data  function  f(x)  for  a  <x  <b  can  be  represented  by 
a  Fourier  series  of  the  form 

o»  _  m  - 

f(x)  «  Aq  +  ^2  Afccos (k  ^x)  +  ^  Bksin(k  |-fx)  (1) 

k-1  k-1 


where  D  -  b  -  a,  provided  that  f(x)  is  finite  for  x  in  (a,b)  and  has  a  finite 
number  of  discontinuities.  In  practice,  the  infinite  series  is  truncated  af¬ 
ter  the  N-th  harmonic  for  a  chosen  N  to  make  the  series  finite.  This  gives 
the  approximation 


f(x) 


a  Ao  +  ^  Afccos  (k  I*1  x) 
k-1 


k-1 


Bksin(k  |^x). 


(2) 


The  approximation  is  improved  if  more  harmonics  are  Included  In  equation  (2), 
that  is,  if  a  larger  N  is  chosen;  however,  this  is  not  entirely  true  for  the 
case  of  a  real  dlscrete^data  function.  For  the  case  of  a  real  discrete-data 
function,  there  is  an  optimum  truncation  for  the  Fourier  series  representa¬ 
tion. 

Consider  a  real  discrete-data  function  f(n)  for  n^  <  n  < n2  where  n,  n^  and 
n2  are  integers.  This  function  can  be  represented  approximately  by  a 
discrete-data  Fourier  aeries  of  finite  harmonics,  namely, 

K  K 

f(n)  -  Ao  +  Akcos(k  |pn)  +  ^  Bksin(k  ^pn)  (3) 

k-1  k-1 

where  N  •  n2  -  nj  and 

l  N 

Ao--^-^  f(n)  (4) 

n-1 


l 


Afc  -  -jj  ^2  f(»)co«(k  |pn) 


Bk-|  2  f(n)«in(k  |Jn) 


Note  that  N  la  also  an  integer  which  is  the  number  of  data  points  representing 
f(n)  in  the  Interval  from  nj.  to  n2*  The  following  interesting  and  Important 
fact  exists: 

For  a  set  of  N  data  points  representing  f(n), 
the  optimum  finite  Fourier  series  consists  of 
the  constant  term,  the  fundamental,  and  all 
harmonica  up  to  and  including  N  . 

7 

Thus,  for  the  optimum  representation,  (3)  becomes 

N/2  ,  N/2 

f(n)  -  Ao  Akcos(k  |Jn)  +  ^  BfcsinCk  |-n).  (7) 

k-1  k-i 

III.  VERIFICATION 

Three  of  the  ways  this  fact  can  be  verified  are  by  sampling  process, 
exponential  form  of  the  discrete  Fourier  series,  and  solution  of  simultaneous 
equations. 

A.  Sampling  Process 

Review  the  sampling  process  by  considering  a  sampled  signal  as 
shown  In  Figure  I.  The  following  quantities  can  be  established: 

Record  length  ■  D 

Fundamental  frequency  F  ■  •jj 

Amber  of  data  points  "  N 

D 

Sampling  period  d  ■  — 

1  N 

Sampling  frequency  fg  ■  ^  ^  ■  NF 

1  NF 

Folding  frequency  fy  •  "jfs  “  j" 


X 


Figure  1  •  Sampled  signal. 


Thus,  the  folding  frequency  is  the  harmonic  j  of  the  Fourier  series.  Figure 
2  shows  the  frequency  spectre  of  a  continuous-data  signal  from  which  the 
discrete-data  signal  is  obtained  by  s sapling,  components  of  the  saapled  signal, 
and  the  saapled  signal.  Rote  that  the  spectrua  of  the  saapled  signal  is  peri¬ 
odic  with  a  period  equal  to  fg,  the  s sapling  frequency. 

For  practical  purposes,  the  high  frequency  end  of  the  saapled  signal  is  at 
fy,  the  folding  frequency.  Thus,  the  entire  practical  spectrua  for  the  saapled 
signal  spans  froa  zero  frequency  to  the  harmonic  H.  There  is  no  useful 

2 

frequency  component  above  the  folding  frequency.  In  fact,  if  the  original 
analog  signal  contains  harmonics  above  fy,  they  are  folded  back  to  the  low 
frequency  region  and  are  added  to  low  frequency  components,  causing  the  ali¬ 
asing  effect.  The  accuracy  of  Fourier  coefficients  for  the  base  band  compo¬ 
nents  can  be  impaired  by  the  aliasing  effect.  The  effect  of  aliasing  Is  less 
if  the  sampling  frequency  is  higher,  implying  a  larger  number  of  data  points 
for  a  glvei  length  of  record. 

The  discrete-data  Fourier  coefficients  are  not  necessarily  equal  to  the 
Fourier  coefficients  of  the  original  analog  signal  due  to  the  picket-fence 
effect.  However,  the  finite  Fourier  series  representation  of  the  signal  may 
agree  well  with  the  original  analog  signal  at  sampling  points. 


0 


-'s 


a)  The  original  signal. 


*S 


b)  Components  of  sampled  signal. 


c)  The  resultant  of  sampled  signal. 


Figure  2 .  Frequency  spectre. 
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B.  Exponential  Fora 


Next,  consider  the  exponential  form  of  the  discrete  Fourier  series. 
It  is  well  knownl  that  the  discrete  Fourier  series  representation  of  a  peri¬ 
odic  sequence  f(n)  of  period  N  needs  only  contain  N  conplcx  exponential  terns 
nanely. 


N-l 


f(n)  -  £  F(k)  e 


4  2*  u 

J  5-“k 


k-0 


(8) 


where  the  F(k)'s  are  Fourier  coefficients.  One  sust  be  careful  to  note  that 

N 

equation  (8)  contains  only  the  constant  tern  and  all  harmonics  up  to  y  inclu¬ 
sive,  but  not  up  to  the  N-th  one.  This  is  due  to  the  fact  that  equation  (8) 
is  equivalent  to 


and 


N/2 

f(n)  -  £ 
k-  - 


F(k) 

N/2 


e 


4  21T  L. 

J  gr-nk 


(9) 


F(k)  -  F(-k)*  -  F(N-k)* 


(10) 


where  denotes  couplex  conjugate.  With  the  help  of  Euler’s  formulas,  equa 
tion  (9)  can  be  reduced  to  the  fom 


f(n) 


N/2 

k-1 


sw 


(U) 


which  does  not  have  hamonlca  above 


N 

T 


^A.V .  Oppenheln  and  R.  W.  Schafer,  Digital  Signal  Processing. 
Prentlc-Hall,  1975,  p.  88. 


C.  Slailtaneous  Equations 


Finally,  consider  aquation  (11)  for  n  •  1  to  N  which  gives  a  sat  of 

N  equations.  On  the  right-hand  side  of  each  of  these  equations  there  are  N 

coefficients  to  be  determined.  They  are  A o,  A^  for  k  ■  1  to  £  and  B^  for  k  - 

N  * 

1  to  y  -1.  Note  that  Bfl/2  ■  0.  Hence  there  are  N  equations  to  be  solved  for 

N  unknowns,  giving  a  set  of  unique  solutions.  In  other  words,  the  sequence 
f(n)  can  be  perfectly  represented  by  a  finite  discrete  Fourier  series  con¬ 
taining  the  constant  tern  and  all  harmonics  up  to  ^  inclusive,  where  N  Is  the 
nuaber  of  data  points.  * 

IV.  RESULTS 

A  numerical  example  Is  given  to  illustrate  the  fact  discussed.  Consider 
a  set  of  36  data  points  representing  the  errors  of  a  resolver  at  36  equally 
spaced  angular  positions  in  the  range  from  0  to  360  degrees,  as  shown  in  Table 
1.  The  coefficients  for  the  corresponding  discrete  Fourier  series  are  shown 
In  Table  2.  Notice  that  Bjj/2  ■  B^g  ■  0,  and  values  of  the  coefficients  are 

N 

symmetrical  about  the  y  ■  18  point.  Table  3  shows  the  Fourier  representation 

RMS  error.  Notice  the  error  is  minimum  when  the  series  contains  only  the  con¬ 
stant  term  and  all  harmonics  up  to  £  inclusive.  In  fact,  this  minimum  value 

2 

should  be  zero.  The  discrepancy  Is  caused  by  the  finite  word-length  effect  of 
the  computer  used.  Figure  3  depicts  the  RMS  error  versus  the  included 
harmonics. 
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TABLE  2.  Fourier  Coefficient*  of  the  Resolver  Error  Cherect eristics. 


HARMONICS 

COS  COEFF 

SIN  COEFF 

EXP  COEFF 

K 

A<K> 

BOO 

COO 

0 

107.383333333 

107.583333333 

1 

-80.6384477067 

28.2839909441 

85.4473711772 

2 

-120.756430444 

-48.9906110651 

138.696332025 

3 

-.332407730167 

-10.5587235027 

10.5639546361 

4 

46.2600631468 

-58.4540750361 

74.3505759837 

5 

4. 3000044443 

-2.03413833927 

4. 76582330616 

6 

3.36111111122 

-3.31221413737 

6.33288164731 

r 

.375601020306 

-4.26831052972 

4.28430473076 

8 

131. 133323603 

-42.6963543421 

137.956363142 

9 

- 1 . 33333333333 

-.333333333333 

1 . 36338998124 

10 

1 . 30414316343 

-6.43221895044 

6 . 6  3  0  4  4  7  0  8  0  7  4 

1  1 

£. 07084910347 

—  6. 7 4962 1 344 7 

7. 06015628335 

12 

-3. 58333333322 

-4. 18378943159 

9.54957823553 

13 

.320306230261 

2.27492947856 

2.4133743756? 

14 

-5.21437939322 

1 . 56652380156 

5.44460712731 

15 

-2.33423391661 

. 892056836222 

2. 97132781856 

1 6 

-2.93343384339 

-14.2644932841 

14.5670472299 

17 

-7.94641400322 

-1 . 87602743122 

8. 16436218549 

18 

3.61111111111 
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3.61 1  1  1  1  1 1  1 1  r 
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-7.94641400322 

1.87602743122 

8'.  1643621834? 
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-2. 95343884389 

14.2644932841 
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21 

-2. 33423391661 

-. 892056336222 

2. 97132731356 

"T* 

-5.21437939322 

- 1 . 56632380 1 56 

5.44460712731 

23 

.320506239261 

-2.27492947356 

2. 4183743756? 

24 

-3.38333333322 

4. 18373945159 

9. 54957323553 

25 

2.07034918547 

6. 74962132478 
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26 

1.30414316345 
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6 . 6»044 1 0 8 O 7  4 

-1.83333333333 
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1 . 36333998124 

23 

1 3 1 . 1 83325698 

42. 6963543421 

137.956363142 

20 

. 37560 1 920586 

4.26331052972 

4. 28430473076 

30 

3.86111111122 

3.31221413757 

6. 83238164731 

31 

4.3099O44445 

2.83413335927 

4. 76582380616 

32 

46 . 269963 1463 

53.4340730361 

74.5505759887 

33 

-.332407750167 

10.3337233027 

10. 5639546061 

34 

-129. 736430444 

48. 99061 10651 

138.696332025 

35 

-88. 6384477067 

-23.2839909441 

35.4473711772 
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TABLE  3.  Fourier  Representation  RMS  Error  Using  Different  Numbers 
of  Harmonic  Terms. 


RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 
RMS  ERROR  USING 


10  HARMONICS*  13.3418581462  ARCSEC 

11  HARMONICS*  14.3068843482  ARCSEC 

12  HARMONICS*  12.8394893114  ARCSEC 

13  HARMONICS*  12.7231019045  ARCSEC 

14  HARMONICS*  12.1287404372  ARCSEC 
13  HARMONICS*  11.9433736037  ARCSEC 

16  HARMONICS*  6.04917499303  ARCSEC 

17  HARMONICS*  1.80333333316  ARCSEC 

18  HARMONICS*  1.30533333388  ARCSEC 

19  HARMONICS*  6 . 049 1 7499673  ARCSEC 

20  HARMONICS*  11.943373609  ARCSEC 

21  HARMONICS*  1 2 . 1 237404603  ARCSEC 

22  HARMONICS*  1 2 . 723 1 0 1 9076  ARCSEC 

23  HARMONICS*  1 2 . 8394893 1 46  ARCSEC 

24  HARMONICS*  14.306334833  ARCSEC 

25  HARMONICS*  15. 3413531503  ARCSEC 


Figure  3.  RMS  error  versus  harmonics  included. 
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